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ABSTRACT 

In a previous paper, we outlined a new Bayesian method for inferring the properties of extended 
gravitational lenses, given data in the form of resolved images. This method holds the most promise 
for optimally extracting information from the observed image, whilst providing reliable uncertainties 
in all parameters. Here, we apply the method to the well studied optical Einstein ring 0047-2808. Our 
results are in broad agreement with previous studies, showing that the density profile of the lensing 
galaxy is aligned within a few degrees of the light profile, and suggesting that the source galaxy (at 
redshift 3.6) is a binary system, although its size is only of order 1-2 kpc. We also find that the mass 
of the elliptical lensing galaxy enclosed by the image is (2.91±0.01)xl0 11 M Q . Our method is able 
to achieve improved resolution for the source reconstructions, although we also find that some of the 
uncertainties are greater than has been found in previous analyses, due to the inclusion of extra pixels 
and a more general lens model. 

Subject headings: gravitational lensing - methods: data analysis - methods: statistical 



1. INTRODUCTION 

Gravitational lenses have long held the promise 
of being natural telescopes, probing structural scales 
smaller than the resolution of current technologies 
while simultaneously revealing the distribution of mat- 
ter within the lensin g galaxy (For reviews of grav- 
itational lensing, see Schn eider. Ehlers. fc Falcol Il99"l 
iKochanek fc Wambseanssl l2004f) . In practice, this 
promise can only be truly realised in systems in which the 
lensed images are extended, as these can offer substan- 
tially more constraints than multiply imaged point-like 
quasars. For this reason, significant attention has been 
paid to extended sources lensed by intervening galaxies 
to pr oduce roughly circular images know n as Einstein 
rings l|Kochanek. Keeton. fc McLeodll200l|) . 

A long-standing issue has focused upon the ques- 
tion of what is the best approach to invert such ex- 
tended gravitationally lensed images, providing the un- 
lensed source brightness distribution and mass in the 
lensing galaxy. Addressing this question has spawned 
a number o f seemingly distinct m ethods, such as the 
ring Cycle l | K ochanek _et alJ fl989f> . Semi-Linear Inver- 
sion llWarren fc Dvd 12 003) and Genetic Algorithms 
l)Brewer fc Lewisf l^OO.'i'l. Given the extended nature of 
the source, reconstructions have been based upon a pixel- 
lated source plane, so as to make minimal assumptions 
about the source brightness distributions. However, due 
to the ill-posed nature of such an inversion, authors have 
tended to use low numbers of pixels, and/or regulariza- 
tion. In a previous contribution, we demonstrated that 
other approaches can be unified in terms of a Bayesian 
interpretation, with each corresponding to the use of par- 
ticular, often u njustified, assumptions about the nature 
of the source (. Brewer fc Lewis) 120061 hereafter BL06). 
Furthermore, BL06 presented a general Bayesian ap- 
proach to the question of gravitational lens inversion, us- 
ing realistic prior distributions and Markov Chain Monte 
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Carlo (MCMC) methods to recover the properties of the 
lensing system. 

In this present contribution, the approach detailed in 
BL06 is applied to the well-studied optical Einstein ring 
0047-2808. In Section the details of this system, and 
the data employed, are presented, while Section|3]details 
the approach to the problem. Section 0] outlines the re- 
sults of this study, with a comparison to other techniques, 
while the conclusions are presented in Section [SJ 

2. ER 0047-2808 

2.1. Background 

The optical Einstein ring 0047-2808 was identified 
serendipitously in a survey of massive elliptical galax- 
ies at z ~ 0.5 via the identification of an anomalous 
emission line, with subsequent imaging revealing a ring 
of high redshift (z = 3.6) emission superimpo sed upon 
a z = 0.49 galaxy fWarren et all 119961 IT999|) . Initial 
estimates suggested that the source was magnified by a 
factor of ~ 17 and hence detailed spectroscopy is able 
to prob e this high redshift. kiloparsec scale star-forming 
galaxy l| Warren et al-l UggS). 

2.2. Observations 

Since its discovery l)Warren et al.lll996D . the optical 
Einstein ring 0047-2808 has been the subject of sev- 
eral different studies whose goal was a full lensing in- 
version and reconstruction of the source brightnes s dis- 
tribution iWarren et all 119991 iKoopmans fc Trel 120031: 
iWavth et al.ll2005t iDve fc Warrenll2005[) . The more re- 
cent of these have focused upon observations obtained 
with the Wide Field Planetary Camera 2 (WFPC2) on- 
board the Hubble Space Telescope and given its superior 
resolution over previous ground-based observations, this 
data is the subject of this current contribution. 

The raw pixel scale of the WFPC2 CCD used in imag- 
ing 0047-2808 is 0.1 arcseconds, but the image used in 
this study is comprised of four interlaced dithered im- 
ages, giving a pix el resolution of 0.05 arcseconds (see 
IWavth et al.l 120051 for details). The foreground lensing 



2 





1.5 




1.0 


ds) 


0.5 


£= 




O 




O 




03 





CO 


o 

: — 




< 


-0.5 








-1.0 




-1.5 



.-■■■L r -;r 



■ - 



.5 -1.0 -0.5 0.5 1.0 1.5 
x (Arc Seconds) 



Fig. 1.— The HST/WFPC2 image of the optical Einstein ring 
0047-2808 employed in this study. The image has a side length of 
3.05 arcseconds and comprises of 61x61 pixels. As noted in the text, 
the ligh t distribut ion o f the foreground galaxy has been subtracted 
fsee lWavth et all2005l) : this foreground light distribution increases 
the noise in the central regions of this image. 



galaxy has been subtracted, after fitting its brightness 
distribution with a Sersic profile. A region of 61x61 pix- 
els, centred on the resulting lensed image and encom- 
passing a region of 3.05x3.05 arcseconds, was extracted 
and employed in this study; this image is presented in 
Figure [I] 

As with lWavth et~aTI (f2005^ , a noise frame covering the 
same region and including the various noise components 
(such as that from the light distribution of the lensing 
galaxy) was also employed. It should be noted that this 
study did not recreate the image dithering procedure that 
was undertaken with the observed system, rather it is 
assumed that the image presented in Figure ^ is a true 
representation of this lensing system at a resolution of 
0.05 arcseconds. To this end, the point spread function 
was generated wi th the TinyT im algorithm to match this 
resolution scale i|KristJ IT995|) . Although in reality the 
noise values between neighbouring pixels are correlated, 
we will use an uncorrelated likelihood function, which is 
computationally much more managable and will lead to 
slightly more conservative results. 



3.1. 



3. APPROACH 
Parameter Estimation 



The basis of our method was presented in BL06. For a 
more general introduction to Bayesian methods, includ- 
ing Mar kov Chain Mon te Carlo techniques, see the text- 
book bv lGregorvl (|200 •">;•. Assuming a particular form for 
the lens model, and background information or assump- 
tions J, we write down the joint probability distribution 
for the lens and source parameters (denoted collectively 
by L and s respectively) given the data D and the prior 
information and assumptions /: 



p(SL\DI) = 



p(SL\I)p(D\SLI) 
P(D\I) 



(1) 



The denominator does not depend on the source and 
lens parameters, so is part of the normalisation constant 



for the joint probability distribution. The first factor 
in the numerator is the joint prior probability density 
of the lens and source parameters, and the second fac- 
tor is the likelihood function (probability density of the 
data that was actually observed, given the parameters). 
Commonly, the source is pixellated into m pixels and 
is described by a vector s of pixel intensities, and the 
observed image is also pixellated (with n pixels), and 
is represented by a vector O. In this case the image 
predicted by a source is calculated by a matrix multi- 
plication: I = Ls. The matrix L depends on the lens 
parameters L and the point spread function. Since it is 
always possible that the noise level values {a{\ are not 
completely reliable, we also introduced a n extra noise pa - 
rameter a, to be estimated from the data l)Gregorvl 2005'l. 
If a good fit is possible with our model, a will be esti- 
mated to be small. If a good fit is not possible, a will be 
estimated to be large, and the uncertainties in all other 
inferred parameters will also be increased. Thus, we can 
use the inferred value of a as a kind of alternative-free 
test of whether our lens model is adequate. Under the 
usual assumptions (a Gaussian probability distribution 
for the error in each pixel of the image, and logical in- 
dependence of the prior information about the lens and 
the source), we have the PDF 



p(sLa\DI)cxp(s\I)p(L\I)p(a\I)p(0\sLaI) 
(xp(s\I)p(L\I)p(a\I) x 



n 



i 



(O; 



(2) 



where the right hand side depends on the lens parame- 
ters through the matrix L in some complicated nonlinear 
fashion. Therefore, there is no hope of finding marginal 
distributions analytically, so we will use a Markov Chain 
Monte Carlo algorithm to generate random samples from 
this posterior distribution. 

3.2. Lens Model 

For its convenient properties, we chose an softened 
power-law ellip tical potential (SPEP) for the lens model 
(Barkana 1998). Elliptical potentials have the advantage 
of faster computation of the deflection angles than for 
elliptical projected mass distributions, and great flexibil- 
ity in the range of possible lenses is achievable with a 
handful of parameters. Hopefully, this flexibility is suf- 
ficient to ensure that we won't make any overconfident 
conclusions caused by the fact that real lenses are not ex- 
actly described by some parametric model. The primary 
reason for choosing this model was computational speed, 
which is an important consideration because our MCMC 
algorithm requires a lot of likelihood evaluations. 

The 2-D lensing potential we used was 



4>{x,y) = -FJ 
1 



where 



R= V~r 



x 2 q ■ 



y 2 /q 



The gradient of <fr gives the deflection angles as 
a x {x,y) = ^ = {bqx)B?- 2 



(3) 
(4) 

(5) 
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a y (x,y) = ^ = (by/q)E?- 2 (6) 

The parameter b is the overall strength of the lens, 
and q is the ratio of the ellipse's major and minor axes. 
q ~ 1 indicates an approximately circularly symmetric 
lens. Furthermore, (x,y) is position in the lens plane, 
in a coordinate system aligned with the lens and centred 
at the centre of the lens. Thus, in fitting the lens to an 
observed image, we have 7 free parameters: b, q, (x c ,y c ) 
(the position of the centre of the lens), 7, a slope param- 
eter, the core radius r c and 9, the angle of orientation 
of the lens mass distribution with respect to the axes 
defined by the image. If q — 1,7 = 1 and r c = 0, this 
model reduces to the common isothermal sphere model, 
for which b is the Einstein Radius. It is not clear whether 
(x c , Vc) were considered known a priori in previous inves- 
tigations; here we will consider them partially known by 
assigning a weakly informative prior probability density. 
Using the scaled lens equation with source and image 
plane positions measured in arc seconds corresponds to 
a choice of the units of mass, as being that which would 
make the Einstein radius equal to 1 arc second. 

It is usually the case that an image can be fitted with a 
variety of differently parametrized lens models. However, 
for most purposes we are more interested in the values 
of the parameters themselves (e.g. how much mass is 
there, how elliptical is the density profile, is it aligned 
with the light profile of the lensing galaxy, etc) and little 
would be gained by comparing different choices for the 
parameterisation of the mass model. However, we sus- 
pect that, if very high resolution images were available, 
simple lens models would be inadequate and we would 
need to use a nonparam etric mass model, similar to that 
used by Marshall (2005) for weak lensing. It is difficult to 
determine at what point this would be neceessary, how- 
ever, if our "extra noise" parameter a was inferred to 
be large then this would be a definite indication that a 
different mass model is required. 

3.3. Lens Parameter Priors 

In order for the problem to have a definite solution, 
we need to introduce prior probability distributions for 
the parameters of the lens model, to describe our prior 
knowledge or ignorance about their values. In principle, 
we could consider the image of the foreground galaxy as 
providing prior information, however we did not carry 
out any sophisticated analysis along these lines, because 
the final results depend only very weakly on the specific 
choice of functions. The prior probability densities we 
used for each parameter are shown in Tabled These are 
only intended to capture vague prior information about 
the range of values that the parameter could plausibly 
take. For the position of the centre of the lens model, 
we assumed this was near the centre of the light profile 
of the foreground galaxy, but the large prior standard 
deviation of 0.2 means that this is a weak assumption. 

3.4. Source Prior 

Since the source may have complex structure, we need 
to introduce enough pixels to capture this in detail. We 
chose to use a 48x48 pixel grid to represent the source. 
The range of the source plane that was covered was 0.6 



TABLE 1 

Prior probability densities for the lens model parameters, 
and also the extra noise parameter a. 



Parameter Prior Distribution 

~b Exponential, mean 

q Normal, mean 1", SD 0.2" 

r c Exponential, mean 0.5" 

7 Normal, mean 1, SD 0.2 

x c Normal, mean 0.092", SD 0.2" 

y c Normal, mean 0.164", SD 0.2" 

9 Uniform, between and 2tt 

ct Jeffreys' uninformative prior ocl/cr 



arc seconds across, so the resolution over the source plane 
was a factor of 4 greater than the resolution of the image. 
If the image strongly constrains part of the source, this 
will be reflected in the posterior samples. This may or 
may not occur (depending on the quality of the data), 
but using large numbers of pixels at least keeps the op- 
tion open. Since we are not using an optimisation based 
approach, overfitting does not occur. 

In BL06, we discussed the shortcomings of regulariza- 
tion based priors for astronomy, and suggested simple 
priors which should be more realistic. In this paper, we 
use a prior which is somewhat more complex than that 
described in BL06. The prior probability distribution 
describing the range of source models which we consider 
plausible was constructed as follows. Since the source 
is mostly blank, rather than trying to infer each source 
pixel value directly, we imagine that the source is gen- 
erated by "monkeys" throwing N "atoms" of intensity 
onto the (initially blank) source plane. Each atom has 
four attributes: An intensity B, a discrete position 
indicating which pixel the atom landed in, and a width 
parameter. The prior probability density of the intensity 
B of an atom is taken as exponential with mean \x (cho- 
sen to match the typical brightness scale of the image), 
and the prior distribution of the discrete position 
of each atom is uniform. T his is similar to the Massive 
Inference prior suggested bv lSkilliriel l)1998f) for positive, 
additive distributions such as surface brightness distribu- 
tions. It has the desirable property that the prior proba- 
bility distribution for the amount of integrated light over 
any macroscopic region R of the source plane is inde- 
pendent of how we choose to subdivide that region into 
pixels. 

Allowing the atoms to have variable size allows us to 
reconstruct plausible astrophysical sources with only a 
small number of atoms, negating the problems caused by 
an entropic prior, which tends to make the reconstructed 
sky too bright, possibly also causing other parameters to 
be inferred incorrectly (BL06). MEM may be suitable 
for photographs, but for astronomy, Masslnf is more ap- 
propriate. Also, there is some prior expectation that the 
source should contain correlations, so that bright pixels 
are more likely to be near other bright pixels. To take this 
into account, we allowed three different types of atoms 
to make up the source - single pixel sized atoms, and two 
"fuzzier" types of atoms that spread light over several 
pixels. The effect of this variable atom size is to reduce 
the impact of our particular choice of pixellation scale, 



4 



by allowing finer structure in those parts of the source 
where it is justified by the data. We chose to discretize 
into three distinct types of atom to reduce computational 
load, because if we had allowed the atom width to be a 
continuous variable it would have taken much more com- 
putation to produce the pixellated source for lensing. 

In our MCMC simulations, we modified the source by 
manipulating the atoms (positions, intensities, widths, 
and number) directly, rather than by modifying source 
pixel values directly. This helps, because an increased ac- 
ceptance rate can be achieved using proposal transitions 
that move an atom from one pixel to its neighbour; since 
this will only have a slight effect on the predicted image, 
the acceptance probability will be fairly high. Also, most 
proposal moves for the Metropolis algorithm will be con- 
centrated around modifying the bright parts of the source 
(where the atoms are), so little time is wasted adjusting 
pixels where the source is mostly blank. 

3.5. MCMC Sampling 

A feature of the lens inversion problem is the fact that 
it is much faster (by a factor of ~ 50) to compute the 
image for a new source, with a fixed lens model, than it is 
to compute the image for a new lens parameter set. This 
is because, to calculate the image after modifying only 
the source, a matrix multiplication (by a sparse matrix) 
is required. However, to calculate the new image after 
a change in the lens parameter values requires that we 
calculate the sparse lensing/blurring matrix L, which re- 
quires a large ray tracing calculation, firing 100 rays per 
pixel in order to get an accurate approximation to the 
lensing matri x. A modifie d Metropolis method has been 
developed bv INeall l)2004[) which is designed to improve 
the sampling when some variables are "slow" and oth- 
ers "fast" , so we implemen ted this in our code. Parallel 
tempering l|Gregorvll2005j) was also used to ensure effi- 
cient exploration of the parameter space. The MCMC 
simulation was run several times with different random 
number seeds, in order to check convergence. The results 
were deemed to be reliable enough if the variance of the 
error bars returned was significantly less than their size, 
which was the case. 

4. RESULTS 
4.1. Source 

Figure [21 shows two results from the MCMC simula- 
tion (i.e. two samples from the posterior distribution). 
The lensed and blurred images are also shown, and the 
residuals between the images and the observed image are 
consistent with pure noise, in terms of \ 2 . If the num- 
ber of degrees of freedom (data points minus effective 
number of free parameters) was well defined, the image 
would not match the observed one to the extent that a 
frcqucntist goodness of fit test would demand. However, 
since we are not using optimization, these common tests 
do not apply. 

Individual reconstructions tend to have their flux con- 
centrated into a small number of pixels (across a physical 
length scale of 1-2 kpc 1 ), but the reliable conclusions are 
the ones which are reproduced across most of the sources 
in the sample. Nine more samples are shown in Figure|3| 

The assumed cosmological parameters were f2 m = 0.27, Q\ = 
0.73 and H = 71 km s" 1 Mpc" 1 . 
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Fig. 2. — Two samples from the posterior distribution of sources, 
along with their lensed blurred images and the normalized residu- 
als between the model and the data. The source plane is 0.6 arc 
seconds across, corresponding to a physical length scale of ~ 7 kpc. 



In Figure 21 the posterior mean source is displayed. 
This may be taken as a final estimate of the source (if 
a single estimate is required), and when this is lensed, 
it is consistent with the image to within the size of the 
noise error bars. The posterior samples and the mean 
source both show strong evidence for complex structure 
in the upper right region of the source, with a compo- 
nent protruding out at right angles to the main compo- 
nent. Other authors have sugges ted that this may be 
a system of intera cting galaxies (jDve fc Warrenl 120051: 
iWavth et al J 12005). This feature is present in the vast 
majority of the source samples, indicating that the evi- 
dence for two components is strong. 

In principle, to decide whether or not this is a double 
source, we should measure the posterior probability that 
the source is double, by counting the fraction of samples 
which are double. But this is problematic - since we are 
using a nonparametric source model, it becomes hard to 
define what is meant by a double source. Despite this 
difficulty, the question can be answered by simply ob- 
serving the sequence of source models and noting that 
the component at right angles to the main light source is 
persistent across the vast majority of the posterior sam- 
ples (Figure 0), indicating a robust conclusion. This is 
also demonstrated by the significance maps in Figured 

In the primary component of the mean source map, 
some substructure is visible, suggesting that we may have 
resolved extra detail in the core of this galaxy. However, 
the low value of the signal to noise ratio (~ 1) casts doubt 
on this conclusion. Individual pixels are highly under- 
constrained in our approach, it is the integral properties 
of the source that are robustly measured. As a check, 
we calculated the posterior probability that both pixels 
(20,27) and (20,29) are brighter than pixel (20,28), and 
found this probability to be 0.48, and hence that this 
question cannot be answered with the current data. Also, 
there is some weak evidence for more patches of the light 
to the left of the primary component, but better data 
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Fig. 3. — More sample sources from the posterior distribution. 
The diversity across the samples is an indication of the uncertainty 
that we have about the nature of the source. 
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Fig. 4. — The mean source reconstruction (top left) serves as a fi- 
nal estimate of the source, with the two components clearly visible. 
The next two panels are the marginal posterior probabilities that 
each pixel is greater than zero, and greater than 500 respectively, 
and so quantify the amount of evidence for each part of the source. 
The bottom right plot is the ratio of the mean flux to the standard 
deviation, showing that we have the most information about the 
brightest parts of the source. 



would be required to firmly establish this. 

In Figure the two possible components of the mean 
source are plotted separately, and the images of them 
are shown. The first component contributes to three of 
the images, but causes a large extended region of light on 
the left of the image. The second component is needed in 
order to explain the additional, more concentrated com- 
ponent of the leftmost image, as well as the presence of 
the fourth image. A si milar conclusion was r eached by 
IWavth et all l)2005j) and lDve fc Warrenl l)2005l) . however, 



our simulations have found evidence for a structure that 
is long and narrow, on the scale of the pixel size that we 
have used. Since other authors have used larger pixels, 
they have found that the best fitting source is a binary 
source in which both components are blobby, whereas 
our results suggest a needle-like structure for the second 
component, with a width of ~ 200 pc. 

4.2. Lens Parameters 

The estimates of the lens model parameters are shown 
in Tabled (the units, where appropriate, are arcseconds, 
except for 8 which is i n radians). These e stimates are 
similar to those found bv lWavth et alJ 1)20051 ) using a sim- 
ilar (PIEP) lens model, but our uncertainties are greater, 
by a factor of ~ 1-4. 

There are several factors contributing to this increased 
uncertainty. One is due to the fact that we have used 
more pixels than others, making weaker a priori assump- 
tions about the source's structure. We have also used 
a different, more reasonable prior for the source, tak- 
ing into account the fact that astronomical sources are 
usually localized bright patc hes with a dark ba ckground. 
Regularization based priors ijSuvu et al.ll200fft may not 
take this into account. 

Another possible cause for the increased uncertainty is 
the approximations that have previously been used, such 
as a Gaussian approximation to the posterior PDF about 
its peak. This is typically a good approximation if the 
data constrains the parameters very well, but this is not 
usually the case with "non-parametric" (actually many- 
parametric) reconstructions. Finally, we did not fix the 
centre of the lens mass distribution at a particular point, 
but included this as an extra pair of parameters to be 
inferred from the image, and also included the 7 and r c 
parameters. 

To infer any properties of this lensing system (with cor- 
responding uncertainties), all that is required is to calcu- 
late the desired property for all of the posterior samples, 
and observe the diversity of the answers. As an example, 
we were able to measure the total magnification of this 
system, defined as the flux in the image divided by the 
flux of the source; this value was found to be 17.9 ± 1.7, 
further demonstrating the value of lensing for investiga- 
tions of high redshift galaxies. 

The lens mass distribution was found to be slightly 
elliptical, with q = 0.932 ± 0.006, and the angle of orien- 
tation of the lens (which can be seen from the orientation 
of the critical curves in Figure is close to that of the 
subtracted foreground galaxy ( see IWavth et al] l2005h 

From these results we can also estimate the total mass 
contained within the image. We defined a circle of radius 
1.2 arc seconds centred in the middle of the image and 
integrated the density (proportional to the Laplacian of 
the lensing potential) within this ring, for each lens pa- 
rameter set. Given our assumptions, the mass contained 
within the ring was found to be (2.91 ±0.01)xl0 11 Mm, in 
agreem ent with the value reported bv lKoopmans fc TrevJ 
( 2003), who found a slightly lower value but also used a 
ring of slightly smaller radius. 

5. CONCLUSIONS 

In this paper, we have presented an analysis of the HST 
image of the Einstein ring 0047-2808 using a Bayesian 
procedure. Most of the results are qualitatively con- 
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Fig. 5. — The two components of the reconstructed source, plotted separately to show how they contribute to the image. On the right, 
the estimated source and lensed (not blurred) image are plotted along with the caustics and critical lines for the best fitting lens model. 



TABLE 2 

Results for lens model parameters. Where appropriate, 

the units are arc seconds, except for 9, which is in 
radians. the values found by lwayth et al.i 120051) for the 
similar piep model are also shown. for the 7 parameter, 
the quoted value from wa yth et al.i 120051) was from a 
softened power-law elliptical mass density model. 



Parameter 


Value 


Wavth ct al. (2005) 


b 


1.177 ± 0.016 C' 1 " 


'T) 1.170 ± 0.004 


9 


0.932 ± 0.006 


0.917 ± 0.004 


r c 


0.09 ± 0.07 


N/A 


7 


1.04 ± 0.06 


1.08 ± 0.03 


X c 


0.107 ± 0.006 


N/A 


yc 


0.176 ± 0.008 


N/A 





2.245 ± 0.009 


2.247 ± 0.01 


IT 


Negligible 


N/A 



sistcnt with previous work, leading us to conclude that 
the various different methods available for reconstructing 
gravitational lens systems all give satisfactory results for 
this system, despite the disagreement about the uncer- 
tainties in the parameter estimates. However, this may 
not a lways be the case with other data sets (see lGregorvl 
120051 for an example). However, our use of smaller pixels 
and a more reasonable prior for the source has allowed us 
to reconstruct the source with a greater resolution than 
has been possible with other methods. As a result, we 
have found that the second component of the source is 



not just a single patch of light, but has a narrow structure 
protruding from it; given structure formation in ACDM 
cosmologies, we should expect such objects in the early 
universe to show significant substructures indicative of 
major mergers and accretion. 
Discovered in continuum radio sources ((Hewitt et alJ 
the number, diversity and observati onal detail of 
Einstein rings has continued to grow (e.g. Cohen et alJ 



l2?M llreufc Koonmans1l2tol ICarTTli et al.ll2003fl . Fur 

thermore, optical E instein rings should be uncovere d in a 
upcoming surveys llMiralda-Escude fc Lehailll992ft . par- 
ticularly the SLACS survey llBolton et alJl2006D. and the 
number of cases is st eadily increasing ijSluse et a.U l2003: 
iCabanac et al.ll2005l) . Hence it is vital that progress is 
made in developing and understanding optimal and re- 
liable methods for analyzing them. Model selection will 
probably play an important role in this area in the future 
l)Suvu eta l. 2006), however we doubt that it will be use- 
ful to compare the evidence for one ad-hoc regularization 
formula against another, when it is known in advance 
that neither of them accurately describe the prior infor- 
mation that is available. Despite these concerns, the level 
of qualitative agreement between the results presented 
in this paper and in others show that these debates may 
only be of minor importance in practice. 
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